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, Abstract. Recent observations indicate that kinetic and magnetic energies are not in equipartition 

in the solar wind turbulence. Rather, magnetic fluctuations are more energetic and have somewhat 
steeper energy spectrum compared to the velocity fluctuations. This leads to the presence of the 
so-called residual energy E r = E v — Ef, in the inertial interval of turbulence. This puzzling effect is 
i— | ' addressed in the present paper in the framework of weak turbulence theory. Using a simple model 

of weakly colliding Alfven waves, we demonstrate that the kinetic-magnetic equipartition indeed 

Ogets broken as a result of nonlinear interaction of Alfven waves. We establish that magnetic energy 
is indeed generated more efficiently as a result of these interactions, which proposes an explanation 
<**] , for the solar wind observations. 
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l^: INTRODUCTION 

^ • Magnetohydrodynamics (MHD) provides possibly the simplest framework for describ- 

ing turbulence observed in the solar wind. When compressibility and the dissipation 
effects may be neglected, MHD equations have exact solutions in the form of Alfven 
waves that propagate (in both directions) along the background magnetic field with the 
same constant Alfven speed. The characteristic feature of those Alfven wave packets 
is that their magnetic and velocity perturbations, b(x,f) and \(x 1 t), are either exactly 
aligned or counteraligned with each other: for wave packets propagating against the di- 
rection of the background field, one has v(x,?) = b(x,?), for wave packets propagating 
in the direction of the field, v(x,f) = — b(x,f). Counterpropagating Alfven wave packets 
interact when they overlap, or collide, with each other, e.g. [1,2]. 

The symmetry between magnetic and velocity perturbations in individual Alfven 
waves may suggest that magnetic and kinetic energies should be in equipartition in 
the developed MHD turbulence. Indeed, phenomenological theories of MHD turbulence 
often do not distinguish between the two spectra, assuming (explicitly or implicitly) 
that they are the same. However, recent observational data indicate that kinetic and 
magnetic fluctuations in the solar wind turbulence are not in exact equipartition at 1AU. 
There is slight excess of magnetic energy and magnetic fluctuations have slightly steeper 
Fourier energy spectrum, e.g., [3, 4, 5, 6, 7, 8]. A possible explanation may be that the 
solar wind turbulence is in a transient state at 1AU [9]. Another possibility is that the 
solar wind contains some remnant structures, magnetic field discontinuities or current 
sheets, advected from the sun and not solely generated by the turbulence [10, 11]. Those 
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FIGURE 1. The field-perpendicular Fourier spectra of kinetic, magnetic, and total energies in numerical 
simulations of strong incompressible MHD turbulence with a strong guide field. The Elsasser fields z + 
and z~ are driven independently at large scales, so that no residual energy is generated by the driving 
mechanism. The numerical resolution is 1024 2 x 256, the Reynolds number is Re w 5600. The details of 
the simulations can be found in [12]. 



explanations are interesting and may be plausible. 

Quite interestingly, recently it has been found that incompressible steady-state MHD 
turbulence exhibits a very similar mismatch between the kinetic and magnetic energies 
to that observed in the solar wind [12]. Moreover, it has been discovered that mismatch 
between magnetic and kinetic energies exists even in the case of weak MHD turbulence, 
that is, turbulence of mostly independent Alfven waves where equipartition between ki- 
netic and magnetic energies is generally expected [13, 14]. These developments suggest 
that breakdown of kinetic-magnetic equipartition may, in fact, be a fundamental prop- 
erty of Alfvenic turbulence (rather than an artifact or some nonuniversal plasma effects), 
which can be understood in the framework of incompressible MHD (see Fig. 1). This 
effect may be complementary to the previously proposed explanations for the solar wind 
observations [9, 10, 11]. 

The difference between kinetic and magnetic energies can be characterized by the 
so-called residual energy, E r = J (v 2 — b 2 )d 3 x = J (z + ■ z~)d 3 x, where z^ 1 = v ± b. Early 
studies of MHD turbulence based on some closure assumptions indeed indicated that 
kinetic and magnetic energies can be different [15, 16, 17]. However, the question 
remains relatively poorly explored since most of the models of MHD turbulence are 
based on the quantities like total energy E = \ f(v 2 + b 2 )d 3 x and cross -helicity H c = 
J(\-b)d 3 x, or their combinations, the so-called Elsasser energies E^ = J(z ± ) 2 d 3 x, 
neither of which contains the information about the residual energy. In what follows we 
present a simple model illustrating how residual energy is generated in a fundamental 
nonlinear process of MHD turbulence - collision of counter-propagating Alfven waves. 



A MODEL FOR RESIDUAL ENERGY GENERATION 



To formulate the model, we use the standard representation of the incompressible MHD 
equations in terms of the Elsasser variables, 



where the Elsasser variables are defined as z = v ± b. In these formulas, v is the 
fluctuating plasma velocity (the equations are written in a frame with zero mean flow 
velocity), b is the fluctuating magnetic field normalized by y/4npo, \ A = Bo/y/4Wpo 
is the Alfven velocity corresponding to the uniform magnetic field Bo (so that the total 
magnetic field is B = Bo + b), P = p/po + b 2 /2 is the total pressure that includes the 
plasma pressure, p, and the magnetic pressure, p is the constant mass density, and we 
neglect driving and dissipation terms, e.g. [2, 18]. 

An important fact can be derived from Eq. (1). If we consider field configurations 
where z T (x, t) = 0, then an arbitrary function z^ 1 (x, t) = F ± (x ± \ A t) is an exact solution 
of the nonlinear equations (1). This solution represents a non-dispersive Alfven wave 
propagating along the direction of the guide field with the velocity T^A- Nonlinear 
interactions are thus the result of interactions between counter-propagating Alfven wave 
packets. If during one interaction (one crossing time) the deformation of the wave 
packets is significant, then the turbulence is called strong; if the deformation is weak, 
the turbulence is called weak. Weak turbulence can be studied using perturbation theory, 
e.g., [19, 20]. In what follows we study how weak interactions of counter-propagating 
Alfven wave packets can generate residual energy, even if the residual energy is absent 
originally. 

Let us assume that the strong guide magnetic field Bo is in z-direction. We consider 
our system in a box that is periodic in x and y directions, with the dimensions 2k x lit, 
and sufficiently long in the z-direction. When the interaction is neglected, we assume that 
the waves are shear- Alfven waves with polarizations normal to the background field, and 
that they have the following simple field-perpendicular structures: 



The functions (2, 3), supplemented by Pq = const that follows from the incompressibility 
condition, are the exact solutions of Eq. (1) without the nonlinear terms. This form 
ensures that without the nonlinear interaction, the residual energy is identically zero: 
Zq • Zq = 0. For simplicity, we also assume that the waves propagating in the opposite 
directions are not correlated with each other, that is, (f + f~) = 0. This happens, for 
example, if the functions /+ and /~ are random functions of their arguments with certain 
correlation lengths , which can be interpreted as the typical size of the interacting wave 
packets in the direction of the guide magnetic field Bo. The angular brackets then mean 
averaging over the z-direction (practically, over scales much larger than /||). 

We will consider the case of weak turbulence, that is, ZqIh <C v A l±, where lj_ = 2k 
is the field-perpendicular size of the interacting wave packets. This condition follows 




(1) 



Zo (x,y,z,t) = f + (z + v A t)cos(y)x, 
z (x,y,z,t) = f~(z-v A t) cos(x)y, 



(2) 
(3) 



from the fact that the linear terms in Eq. (1), estimated as vaZq must be larger than 

the nonlinear ones, ZqZq/1±. In addition we assume that L ^> which allows us to 
retain only the field-perpendicular gradient of the pressure term, and ensures that only 
shear- Alfven waves are generated by nonlinear interactions, see, e.g., [20, 21, 22]. Let 
us now assume that interaction is switched on at time t = 0. The solution can then be 
found perturbatively in small parameter £ = z ± /va, that is, we can write: 



+ zf + zf + . 



(4) 



We illustrate the main steps of the derivation using the field z + . It is convenient to 
introduce the new variables t, = z + v^t and % = t; then the equation for z\ takes the 
form: 



d r z+ = -(z -V) z+-V ± Pl 



(5) 



The pressure Pi should ensure incompressibility of z\\ it can be found by taking the 
divergence of both sides of (5). One finds this way that 



= \f + f sin(x)sin(y). 



We can now integrate Eq. (5) with respect to time, keeping t, = const, and get: 
z+ = G+(z,t) 



^ cos (x) sin (y ) x — ^ sin (x) co s (y ) y 



(6) 



(7) 



where 



G + (z,t) = j f + ^)f-^-2v A x)dx. 



(8) 



In a similar fashion, we can derive the first-order correction Zj , which is given by 



H=G-(z,t) 



— ^ cos(x) sin(y)x + ^ sin(jc) cos(y)y 



G-(z,t) = [f + (ri+2v a T)f-(ri)dT, 



(9) 



(10) 



where r\ = z — VAt, x = t, and the integral in (10) is done assuming r\ = const. We can 
proceed this way to the second and higher orders, calculating pressure corrections at 
each step. The expressions for the second-order corrections z^ and Zj are bulky; they 
will not be needed for our current discussion and will not be presented here. 

We can now study the evolution of residual energy after the interaction is switched 
on. Substituting the expansion (4) into the expression for the residual energy, we find to 
the first nonvanishing order 



Ey = (Z + • Z ) = (Z+ • Zj ) = 



(ID 



where the angular brackets mean averaging in the x, y, and z-directions. It can be checked 
that the following terms vanish: (zjj" • Zq ) = (zf ■ zj) = (z* • z^) =0. The averaging in 
formula (11) can be easily performed since the functions /+ and /~ axe not correlated 
with each other. We therefore derive: 

t t 
(G+G-) = J (f+^)f + (i 1 +2v A x))dx J {f-{r])f-^-2v A x))dx. (12) 



We now assume that the distributions of /+ and /~ are uniform along the z-direction, 
that is, their correlation functions depend only on the difference of their arguments: 

(f ± (zi)f ± (z 2 )) = K ± (z l -z 2 ). (13) 

The correlation functions are even, and they decline fast when their arguments 
exceed the correlation length Zy. We also consider Fourier transforms of /^(z), K" ± (z): 

-\-oo -\-oo 

f±iz) = 2^ / ^(QeMitydk, ^(z) = ^ J K- ± (fc)exp(fc)^, (14) 

— oo — oo 

so that their Fourier components are related as: 

±(\f ± (k)\ 2 ) = K ± (k), (15) 

where L is the system size in the z-direction (L 3> /||). Then the integrals in (12) can be 
easily done: 

t t t 

J(f + ($)f + (ri+2v A x)}dx = j K + ^-T]-2v A x)dx = J K + {2v A t-2v A x)dx, (16) 



t t t 

j {f-{^-2v A x)f-{7]))dx = J K-{^-7]-2v A x)dx = J K-(2v A t-2v A x)dx. (17) 



Now, we are interested in times that are much larger that the time necessary for each 
wave to shift by a correlation length, that is, we assume t >> l\\/v A . In other words, we 
are interested in times exeeding a single collision time of two counter-propagating wave 
packets. Then, the integration limits in (16) and (17) can be formally extended to infinity: 

t oo -\-oo 

J K ± (2v A t-2v A x)dx = ^-J K ± (z)dz=^jK ± (z)dz=^K ± (0). (18) 



According to (15), this expression is positive. We therefore obtain that the residual 
energy is negative, and it is expressed through the wave correlation functions as: 

E r = (z+-z-) = -_Ljc+(0)ir(0) < 0. (19) 
This is the main result of our derivation. 



CONCLUSIONS 



An important observation is that the residual energy is always negative, which means 
that the magnetic energy always exceeds the kinetic one. This result, obtained in the 
framework of incompressible MHD, provides a qualitative explanation for the solar 
wind observations. More general discussion indicates that this conclusion is quite robust, 
in that it is independent of the particular form of the initial z + and z~ modes [14]. In 
addition, direct numerical simulations indicate that it holds for strong turbulence as well. 
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